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LONG  TERM  GOALS 

To  develop  a  passive  modality  of  acoustic  thermometry  using  Cross-correlation  processing  of  deep 
water  ambient  noise. 

OBJECTIVE 

Our  previous  research  effort  has  demonstrated  that  coherent  processing  of  shipping  noise 
(100Hz<f<kHz)  in  shallow  water  can  provide  a  totally  passive  means  for  sensing  the  ocean 
environment;  but  optimizing  this  correlation  process  for  this  relatively  high  frequency  band  remains 
challenging  in  a  fluctuating  ocean  enviromnent  with  significant  multipath.  We  hypothesize  that  using 
instead  significantly  lower  frequency  (f<50Hz)  ocean  noise  recorded  in  the  SOFAR  channel  may 
enhance  the  performance  of  this  correlation  process  due  to  the  potential  long-range  spatial  correlations 
(~1  OOkms)  of  this  deep  water  noise  field  and  the  relative  temporal  stability  of  the  SOFAR  channel. 

The  main  objective  of  this  year  research  was  to  assess  the  feasibility  of  passive  (i.e.  noise-based) 
acoustic  thermometry  using  low-frequency  deep  water  ambient  noise  recorded  on  the  global  IMS- 
CTBTO  network  Results  from  this  study  are  expected  to  help  guide  the  design  of  future  cabled  or 
autonomous  ocean  observatory  systems  for  passive  ocean  monitoring. 

WORK  COMPLETED 

Deep  oceans  play  a  major  role  in  absorbing  atmospheric  heat,  thus  measuring  deep  oceans  temperature 
variations  is  necessary  among  others  to  quantify  air-sea  heat  exchanges-  for  instance  to  assess  global 
warming  trends  ( 1 )  and  calibrate  climate  change  models  (2,3).  However,  contrary  to  ocean  surface 
temperatures,  deep  ocean  temperatures  cannot  be  readily  inferred  from  satellite-based  remote  sensing 
methods  ( 4 )  and  are  thus  most  commonly  obtained  in-situ  using  free-drifting  profiling  oceanographic 
floats  but  with  limited  spatial  and  temporal  resolution  as  these  floats  only  cover  the  globe  sparsely  (5). 
Acoustic  thermometry  provides  an  alternative  remote  sensing  modality  for  detecting  fine  variations  in 
deep  ocean  temperatures  based  on  acoustic  propagation  travel  times  between  sources  and  receivers 
inserted  in  the  ocean  ( 6 ).  Furthermore,  a  promising  alternative  to  using  active  acoustic  sources  is  based 
on  noise  correlation  processing,  which  has  successfully  been  used  for  continuously  monitoring,  with 
unprecedented  temporal  resolution,  seismically-active  systems  such  as  fault  zones  (6a)  and  volcanic 
areas  (6b).  More  generally,  experimental  and  theoretical  studies  have  demonstrated  that  an  estimate  of 


1 


the  arrival  time  structure  between  two  receivers  for  a  given  acoustic  or  elastic  medium  (as  detennined 
by  the  Green’s  function  between  the  receivers  locations)  can  be  obtained  from  the  time-averaged  cross¬ 
correlation  of  the  ambient  noise  recorded  by  the  two  receivers.  This  has  been  shown  in  the  fields  of 
ultrasonics  (9,10,10b),  helioseismology  (11,12),  seismology  (6,6a-6f),  structural  health  monitoring 
(13a-13d),  and  elastography  (14).  In  the  context  of  ocean  acoustics,  previous  studies  have 
demonstrated  that  the  integration  time  required  by  this  noise  correlation  method  to  reliably  estimate  the 
discrete  acoustic  paths  traditionally  used  to  infer  temperature  variations  (6)  is  typically  too  long 
compared  to  the  timescale  of  ocean  fluctuations  (15-17).  However,  for  deep  ocean  propagation,  there  is 
an  unusually  stable  acoustic  feature-  the  last  and  most  energetic  arrival  corresponding  to  sound 
propagating  horizontally  along  the  Sound  Fixing  and  Ranging  (SOFAR)  channel-  which  can  be 
estimated  over  a  finite  correlation  time.  The  SOFAR  channel  is  centered  on  the  minimum  value  of  the 
sound  speed  over  the  ocean  depth  and  thus  acts  as  a  natural  waveguide  for  long-range  sound 
propagation  (18).  In  this  study,  we  demonstrate  the  feasibility  of  a  totally  passive  (i.e.  without  using 
any  controlled  active  sources)  method  for  acoustic  thermometry  of  the  deep  oceans  using  only  records 
of  low-frequency  (fi-TOHz)  ambient  noise  propagating  along  the  SOFAR  channel.  Estimates  of  deep 
ocean  temperature  variations  obtained  from  this  passive  thermometry  technique  over  several  years  are 
found  to  be  consistent  with  oceanographic  measurements  obtained  with  classical  free-drifting  floats. 

Here  we  used  continuous  recordings  of  ocean  noise  from  two  existing  hydroacoustic  stations  of  the 
International  Monitoring  System,  operated  by  the  Comprehensive  Nuclear-Test-Ban  Treaty 
Organization,  located  respectively  next  to  Ascension  and  Wake  islands  (see  Fig.  1 A-B).  Each 
hydroacoustic  station  is  composed  of  two  triangular-shaped  horizontal  arrays,  separated  by  L~  130km, 
which  are  referred  to  hereafter  as  the  north  and  south  triads  (see  Fig.  1C).  The  sides  of  each  triad  are 
~2  km  long  and  the  three  hydrophones  are  located  within  the  SOFAR  channel  at  depth  ~1 000m  (18b). 
Examples  of  one-year-averaged  noise  cross-correlations  in  the  frequency  band  [lHz-40Hz]  between 
pairwise  combinations  of  hydrophones  from  the  north  and  south  triads  appear  in  Fig.  1C-D.  From  year 
to  year,  the  nearly  identical  coherent  arrivals  at  time  delays  -corresponding  to  propagation  between 
hydrophone  pairs  of  the  respective  triads  along  the  SOFAR  channel  axis-  build  up  from  noise  sources 
whose  paths  intersect  the  sensor  pairs.  Based  on  the  time-lag  convention  used  for  computing  the 
correlation  waveforms  (18b),  the  coherent  arrivals  shown  in  Fig.  1C-D  occurring  at  positive  time 
delays  are  generated  by  coherent  noise  traveling  northward,  i.e.  from  the  south  triad  to  the  north  triad 
(see  Fig.  IB).  Conversely,  negative  time-delay  arrivals  correspond  to  coherent  noise  traveling 
southward  in  the  opposite  direction  from  the  north  triad  to  the  south  triad  for  the  Wake  Island  site 
(18b).  For  the  Ascension  Island  site,  the  African  continent  blocks  most  of  the  ocean  noise  propagating 
southward  (see  Fig.  1  A)  along  the  line  connecting  the  north  and  south  triads,  so  no  coherent  arrivals 
for  negative  time-lags  could  be  extracted  from  the  noise  correlation  waveforms  at  this  site. 
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Fig.  1.  (A)  Locations  of  the  two  hydroacoustic  stations  (red  dots)  near  Ascension  and  Wake  Islands. 
Map  data  from  Google  Maps.  (B)  Zoomed-in  schematic  of  the  hydrophone  array  configurations  for 
the  Ascension  and  Wake  Island  sites,  which  both  have  a  similar  layout.  Each  hydroacoustic  station 
consists  of  a  northern  and  southern  triangle  array  of  three  hydrophones  (or  triad),  with  each 
triangle  side  equal  to  approximately  2  km.  The  distance  L  between  triad  centers  is  equal  to  126  km 
and  132  km  for  the  Ascension  and  Wake  Island  stations,  respectively.  (C)  Noise  cross-correlation 
waveforms,  averaged  over  three  different  years,  obtained  between  all  9  pairwise  combinations  of  the 
elements  of  the  north  and  south  triads  for  the  Wake  Island.  (D)  Same  as  C,  but  for  Ascension  Island. 
The  beams  shown  by  dashed  lines  in  A  are  centered  on  the  lines  joining  the  centers  of  the  south  and 
north  triads  of  each  hydroacoustic  station  (yellow  line)  and  which  intersect  the  Polar  regions  where 
potential  ice  noise  sources  contributing  to  the  coherent  arrivals  shown  in  C-D  are  located  (18b). 

In  the  low-frequency  band  used  in  this  study  (l-40Hz),  where  most  of  the  acoustic  energy  is  centered 
around  10  Hz,  ambient  noise  propagating  along  the  SOFAR  channel  mainly  originates  from  seismic 
activity  (e.g.  along  major  undersea  fault  lines)  or  ice-breaking  noise  in  the  Polar  Regions  (19-22).  Ice¬ 
generated  ambient  noise  near  the  ocean  surface  in  the  Polar  Regions  (e.g.  due  to  loud  iceberg  cracking 
events  with  levels  up  to  245  dB  re  1  pPa  at  1  m)  can  efficiently  couple  directly  to  the  SOFAR  channel 
due  to  the  natural  change  in  elevation  of  the  minimum  of  the  sound-speed  profile  from  the  sea  surface 
in  the  cryosphere  towards  greater  depths  at  lower  latitudes  where  the  Ascension  and  Wake  Island  sites 
are  located  (21-23).  When  using  averaging  durations  shorter  than  1  year,  the  observed  seasonal 
variability  of  the  amplitude  of  the  coherent  SOFAR  arrivals-  as  show  on  Fig  1C-D-  confirms  the  polar 
origin  of  the  dominant  component  of  the  ambient  noise  generating  these  arrivals.  The  angular  beams 
shown  in  Fig.  1 A  were  extended  all  the  way  to  the  Arctic  and  Antarctic  using  geodesic  paths  to  obtain 
a  simple  estimate  of  the  geographical  area  from  where  ice-generated  ambient  noise  is  likely  to  emanate 
for  each  site.  The  bathymetry  along  the  axis  of  this  angular  beam  indicates  that  the  Ascension  Island 
site  has  a  direct  propagation  line-of-sight  with  the  Antarctica  coastline,  while  at  Wake  Island 
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seamounts  and  islands  partially  block  noise  propagating  from  either  the  Arctic  or  Antarctica.  This 
likely  explains  why  the  amplitudes  of  the  arrivals  exhibit  a  weaker  peak-to-variance  ratio  for  the  Wake 
Island  site  when  compared  to  the  Ascension  Island  site  (compare  Fig  1C  to  ID).  An  additional  factor 
driving  this  difference  in  the  peak-to-variance  ratios  is  that  the  concentration  of  icebergs  -which 
directly  influences  the  ice  noise  source  density  and  therefore  the  amount  of  coherent  noise  traveling 
along  the  SOFAR  channel  up  to  the  hydroacoustic  stations-  is  typically  higher  in  the  South  Atlantic 
than  in  the  South  Pacific  (24). 

Acoustic  thermometry  estimates  average  oceanographic  parameters,  such  as  temperature  variations, 
over  large  scales  by  leveraging  the  nearly  linear  dependency  between  sound  speed  in  water  and 
temperature  ( 25,6)  .  Here  a  totally  passive  version  of  acoustic  thermometry  is  used  that  monitors  the 
fine  temporal  variations  of  the  coherent  SOFAR  arrivals  overall  several  years.  These  arrivals  are 
obtained  using  a  short  averaging  duration  of  1  week,  and  are  used  to  infer  ocean  sound  speed  variations 
within  the  SOFAR  channel  between  the  North  and  South  triads  at  the  Wake  and  Ascension  Island  sites. 
Here  the  SOFAR  channel  extends  approximately  from  390  m  to  1350  m  deep  at  the  Ascension  Island 
site  and  460  m  to  1600  m  deep  at  the  Wake  Island  site,  as  determined  from  the  local  sound  speed 
profiles  and  the  center  frequency  (~10Hz)  of  the  SOFAR  arrivals  shown  on  Fig.  1C-D.  This  short  one- 
week  time-averaging  was  accomplished  with  a  data  analysis  procedure  based  on  using  data  averaged 
over  1  year,  indicated  by  a  time-gated  box  on  Fig  1C-D,  as  a  reference  in  a  signal  processing 
methodology  (i.e.,  beamfonning)  that  directionally  combines  the  coherent  arrivals  collected  between 
triads  (18b).  This  process  allows  us  to  continuously  track  fine  variations  of  the  arrival-  time  structure 
with  a  one  week  temporal  resolution  over  three  years  for  Wake  Island  (Fig.  2A-B)  and  seven  years  for 
Ascension  Island  (Fig.  2C).  These  fine  variations  are  indicative  of  mesoscale  oceanic  fluctuations 
occurring  within  the  effective  depth  extent  of  the  SOFAR  channel  over  the  range  between  triads.  The  1 
week-averaging  duration  was  set  conservatively  to  ensure  a  predicted  error  in  arrival-time 
measurements  (18b)  of  ~4  ms  or  less  at  the  Wake  Island  site  and  ~1  ms  or  less  at  the  Ascension  Island 
site  over  the  whole  observation  period.  The  observed  fluctuations  of  the  arrival-time  errors  across  the 
years  (indicated  by  the  colorbars  on  Fig.  2),  are  primarily  associated  with  Polar  seasonal  variations 
(18b). 
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Fig.  2.  (A)  Temporal  variations  over  3  years  of  the  coherent  SOFAR  arrivals  at  the  Wake  Island  site 
extracted  from  ambient  noise  correlations  using  a  one-week  moving  average.  These  arrivals 
primarily  result  from  ice-noise  originating  propagating  northward  towards  the  Wake  Island  site 
successively  from  the  South  triad  to  the  North  triad.  The  color  associated  with  each  data  point 
indicates  the  estimated  measurement  error  a  in  milliseconds  (18b).  (B)  Same  as  A  but  using  instead 
noise  propagating  southward  towards  the  Wake  Island  site  in  the  reciprocal  direction 
i.e.  successively  from  the  North  triad  to  the  South  triad.  (C)  Same  as  A  but  for  noise 
propagating  northward  towards  the  Ascension  Island  site  successively  from  the  South  triad 
to  the  North  triad.  Due  to  data  availability  a  longer  observation  period  of  ~7 years  was  used 

for  the  Ascension  Island  site. 

Comparing  Fig.  2A  and  B  demonstrates  the  similarity  of  the  arrival-time  variations  obtained  for  both 
the  positive  and  negative  time-delay  arrivals,  which  result  from  different  noise  events  propagating 
northward  or  southward,  respectively,  along  the  same  paths  linking  the  south  and  north  triads  for  the 
Wake  Island  site.  This  symmetry  of  the  arrival-time  variations  (as  confirmed  in  Fig.  S2)  demonstrates 
that  they  are  due  to  reciprocal  changes  in  the  environment ,  such  as  ocean  sound  speed  fluctuations 
induced  by  temperature  changes,  rather  than  non-reciprocal  changes,  such  as  currents,  clock  drift,  or 
other  signal -processing  artifacts.  (26,27).  The  symmetric  arrival  time  variations  could  not  be 
established  for  the  Ascension  Island  site  due  to  the  African  continent  blockage  of  ocean  noise 
propagating  southward.  Instead  a  spectral  analysis  of  the  arrival-time  variations  measured  at  Ascension 
Island  from  northward  propagating  noise, ,  confirmed  a  dominant  period  of  365  days,  consistent  with 
expected  seasonal  variations  of  the  ocean  sound  speed  at  this  site  (18b). 

Finally,  the  fine  arrival-time  variations  with  associated  error  bars  were  used  to  infer  ocean  temperature 
variations  (18b)  averaged  over  the  effective  depth  extent  of  the  SOFAR  channel  over  several  years 
between  the  North  and  South  triads  of  the  Wake  Island  and  Ascension  Island  sites  (see  Fig.  3).  In  Fig. 

3,  overlaid  dots  correspond  to  estimated  monthly  temperature  changes,  and  corresponding  error  bars, 
independently  obtained  from  the  Argo  floats  measurements  (5)  in  the  aforementioned  effective  depth 
of  the  SOFAR  channel  between  the  north  and  south  triads  of  each  hydroacoustic  stations  over  the  same 
time  period.  The  magnitude  of  the  temperature  variations  in  the  SOFAR  channel  estimated  from 
passive  thermometry  and  Argo  data  are  in  good  agreement  over  the  whole  observation  period,  and  the 
measurement  error  bars  from  passive  thermometry  are  significantly  smaller.  An  important 
methodology  difference  between  these  two  estimates  is  that,  given  the  sparseness  of  the  global 
coverage  by  the  Argo  floats,  the  coarse  Argo  measurements  are  spatially  and  temporally  interpolated  to 
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provide  monthly  estimates  of  global  temperature  variations  (5).  On  the  other  hand,  passive 
thermometry  measurements  are  not  interpolated  as  they  result  from  a  true  sampling  of  the  ocean  by 
propagating  sound  and  thus  integrates  temperature  variations  over  the  whole  propagation  path  (here 
-1 30km  between  the  North  and  South  triads  of  each  site)  and  averaging  duration  (here  1  week).  Due 
to  limited  hydroacoustic  data  availability  for  the  newly  constructed  Wake  Island  station,  the 
observation  period  only  covers  3  years,  with  one  data  point  per  month,  and  thus  does  not  provide 
sufficient  statistics  for  the  Argo  data  to  infer  any  significant  long-tenn  trend  in  ocean  temperature 
variations.  On  the  other  hand,  for  the  Ascension  site,  the  SOFAR  channel  temperature  variations 
measured  over  a  longer  duration  of  7  years  from  passive  thermometry  and  Argo  are  found  to  be 
significantly  correlated  with  a  0.78  correlation  coefficient,  indicating  that  this  passive  thermometry 
method  is  sufficiently  accurate  to  detect  mesoscale  variations  of  ocean  temperatures  occurring  within 
the  effective  depth  of  the  SOFAR  channel.  Furthermore,  Fig.  3(B)  indicates  a  warming  of  the  SOFAR 
channel  in  the  Ascension  area  as  inferred  from  the  similar  upward  trend  of  both  passive  thermometry 
(0.017  °C  /year)  and  Argo  (0.013  °C/  year)  temperature  variation  estimates  ( 18b);  this  finding  is 
consistent  with  the  recent  report  of  increased  heat  transport  to  the  deep  southern  Atlantic  ocean  over 
the  same  time  period  (1). 
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Fig.  3.  (A)  Comparison  of  the  deep  ocean  temperature  variations  at  the  Wake  Island  site  estimated 
from  passive  thermometry  using  the  SOFAR  arrival-time  variations  shown  in  Fig.  2A  (blue  line) 
with  free-drifting  profiling  oceanographic  Argo  float  measurements  (grey  dots),  along  with 
corresponding  error  bars  (18b)  (B)  same  as  A,  but  for  the  Ascension  Island  site  using  the  SOFAR 

arrival-time  variations  shown  in  Fig  2C. 

Overall,  the  practical  feasibility  of  passive  thermometry  is  contingent  upon  extracting  coherent  arrivals 
with  sufficient  amplitude  from  noise  correlations  computed  over  the  shortest  possible  averaging  time  T 
in  order  to  accurately  track  temperature  variations  in  a  rapidly  fluctuating  ocean.  Results  from  this 
study  not  only  indicate  that  such  signals  can  be  obtained,  but  also  that  they  can  potentially  be  used  to 
estimate  the  requirements  of  passive  thermometry  (in  terms  of  achievable  measurement  error  for  a 
given  averaging  duration  T)  over  longer  receiver  separation  ranges,  L,  for  temperature  monitoring  of 
deep  oceans  at  basin  scales.  Indeed,  previous  studies  have  established  that  the  N-M  pairs  of  coherent 
arrivals  extracted  between  the  two  hydrophone  arrays  each  containing  a  different  number  of  elements 
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N  and  M,  respectively,  can  be  jointly  processed  (via  the  beamforming  procedure  used  in  this  study)  to 
provide  arrival-time  variations  from  passive  thermometry  with  an  associated  measurement  error 

inversely  proportional  to  V (/ y  26,  28,  29)  , where  L  is  the  separation  distance  between  the 
centers  of  the  two  receiver  arrays.  Note  that  this  simple  scaling  argument  assumes  that  the  statistics 
and  generation  mechanism  of  the  Polar  noise  sources  driving  the  longer  range  results  remain 
comparable  to  the  ones  used  in  this  study.  Therefore,  global  monitoring  of  deep  ocean  temperature 
variations  could  be  enabled  with  high  sensitivity  using  this  passive  thermometry  methodology  by 
leveraging  all  pairwise  combinations  of  a  network  of  hydrophone  arrays  strategically  deployed  around 
the  world. 

CONCLUSIONS 

The  results  of  this  research  effort  confirm  the  possibility  of  passive  acoustic  thermometry  in  the 
SOFAR  channel.  Further  studies  are  required  to  fully  assess  the  potential  of  noise-based  ocean 
monitoring  and  tomography  over  very  long  observation  periods  and  larger  array  separation  distances. 

IMPACT 

It  is  conjectured  that  the  results  of  this  study  could  help  develop  a  totally  passive  means  for  monitoring 
the  ocean  environment  using  only  ambient  noise.  A  potential  scenario  benefiting  from  the  proposed 
methodology  might  include  long-term  deployment  of  ocean  sensing  systems  requiring  minimum  power 
consumption,  and  for  deployments  where  active  sources  are  limited  by  environmental  regulations. 
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